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IS  ABSTRACT 

A  brief  review  of  the  literature  on  blast  wave  theory  is  conducted.  Special 
attention  is  devote-1  to  the  variable  energy  blast  wave  because  of  its  possible 
applicability  to  two-phase  detonation  phenomena  and  the  ignition  of  liquid  mono- 
propcllants  by  shock  waves. 

The  zero  order  solution  for  a  variable  energy  spherical  blast  wave  in  which  the 
total  energy  release  E  is  deposited  proportionally  to  t®,  where  t  is  the  time,  is 
shown  to  be 
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where  R  is  the  shock  radius  and  t^  is  the  duration  of  energy  input.  The  factor  JG 
is  related  to  the  flow  field  via  8 . 

Sample  calculations  of  this  equation  on  the  basis  of  possible  energy  release 
from  ethyl  nitrate  droplets  and  a  comparison  with  experimental  results  are  made. 

It  appears  that  values  of  ft  =  5  and  t^  =  10  usee  give  reasonable  agreemenr.  However 
more  accurate  theoretical  calculations  need  be  made  and  more  experimental  data 
need  be  obtained  before  firm  conclusions  could  be  drawn. 
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NOMENCLATURE 


aQ  =  Speed  of  sound  in  the  undisturbed  gas 

E  =  total  energy  (a=2),  energy/unit  length  (a  =  l),  energy/ 
unit  area  (a=o) 

Ea  =  see  equation  11 

f  =  nondimensional  velocity 

g  =  nondimensional  pressure 

h  =  nondimensional  density 

J  =  constant  defined  in  equation  26 

JQ  =  value  of  J  when  y  ->o 

M  =  Mach  number 

p  -  pressure 

r  =  distance  from  site  of  energy  source 
R  =  shock  radius 

RQ  =  explosion  scale  radius  (equation  34) 

R  =  nondimensional  shock  radius  (equation  42) 
t  =  time  after  start  of  energy  input 
t^  =  duration  of  energy  input 

u  -  velocity 

U  =  shock  velocity 

W  =  proportionality  constant 

=  r/R 


x 


V»* 


y  =  1  /M2 

a  =  geometry  index  =  2  (spherical),  1  (cylindrical),  0  (planar) 
Y  =  ratio  of  specific  heats 
X  =  decay  coefficient  (equation  17) 

Xq  =  decay  coefficient  when  y  ->  o 

p  =  density 

Vlll 


ABSTRACT 


A  brief  review  of  the  literature  on  biast  wave  theory  is  con¬ 
ducted.  Special  attention  is  devoted  to  the  variable  energy  blast 
wave  because  of  its  possible  applicability  to  two -phase  detonation 
phenomena  and  the  ignition  of  liquid  monopropellants  by  shock 
waves. 

The  zero  order  solution  for  a  variable  energy  spherical  blast 
wave  in  which  the  total  energy  release  E  is  deposited  proportionally 
to  t&  ,  where  t  is  the  time,  is  shown  to  be 


where  R  is  the  shock  radius  and  t^  is  the  duration  of  energy  input. 
The  factor  JQ  is  related  to  the  flow  field  via  8  . 

Sample  calculations  of  this  equation  on  the  basis  of  possible 
energy  release  from  ethyl  nitrate  droplets  and  a  comparison  with 
experimental  results  are  made.  It  appears  that  values  of  8  =  5 
and  t^  =  10  (j  sec  give  reasonable  agreement.  However,  more 
accurate  theoretical  calculations  need  be  made  and  more  experi¬ 
mental  data  need  be  obtained  before  firm  conclusions  could  be 
drawn. 
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INTRODUCTION 


The  classical  blast  wave  theory  refers  generally  to  the  propaga¬ 
tion  of  shock  waves  in  a  gaseous  medium  due  to  an  instantaneous 
energy  input  in  an  infinitesimally  small  region  of  that  medium.  The 
energy  input  can  be  in  a  point,  along  a  line,  or  in  a  plane,  with  the 
resultant  wave  and  attendant  flow  fields  being  spherical,  cylindrical, 
or  planar,  respectively.  Interest  in  the  theory  started  in  the  early 
forties  with  the  pioneering  analysis  of  G.  I.  Taylor  (Ref  1), which 
was  not  published  untij  1950.  Since  then  a  very  large  amount  of 
work  on  the  subject  h  .s  been  published  in  the  international  literature. 
For  comprehensive  eviews,  the  reader  is  referred  to  the  works  of 
Sakurai  (Ref  2),  Serov  (Ref  3),  and  Lee,  Knystautas  and  Bach  (Ref  4). 
The  proliferation  c.f  analyses  naturally  has  resulted  in  different  and 
sometimes  confusing  sets  of  nomenclatures.  A  comparison  of  the 
symbols  used  by  important  contributors  for  the  pertinent  parameters 
is  presented  by  Oppenheim  et  al  (Ref  5). 

Early  treatment  of  the  subject  was  based  on  the  assumptions  that 
the  energy  input  is  instantaneous  and  that  transport  properties  are 
unimportant.  With  these  assumptions,  dimensional  analysis  (Ref  3) 
yields  functional  relationships  between  the  energy  input  and  distance, 
arrival  time,  and  velocity  of  the  shock  front  without  complete 
solution  of  the  relevant  conservation  equations.  However,  for  detailed 
characteristics  of  the  flow  field  (i.e. ,  pressure,  density,  temperature, 
and  velocity  of  the  gas  behind  the  shock),  numerical  solutions  are 
generally  required  (Ref  1,2).  Solutions  for  the  flow  field  near  the 
energy  source  region,  where  transport  properties  are  important, 
have  appeared  in  the  literature  (see,  e.g. ,  Bowen  (Ref  6)  where  the 
method  of  asymptotic  matching  technique  is  used).  Problems  in 
which  the  effect  of  the  energy  source  mass  is  important  (piston 
problem)  also  have  been  analyzed  (Goldsworthy  (Ref  7),  Grigorian  et  al 
(Ref  8))„  Also,  Laumback  and  Probstein  (Ref  9)  treat  the  case  of 
point  explosion  in  a  variable  density  atmosphere. 

While  the  assumption  of  instantaneous  energy  input  is  adequate 
for  most  cases  and  certainly  for  time  periods  much  greater  than  the 
input  times,  there  are  processes  in  which  the  energy  input,  though 
very  fast,  can  be  considered  time  dependent.  Examples  are 
discharges,  exploding  wire  phenomena,  accelerating  hypersonic 
body,  and  chemical  energy  release.  An  example  of  the  latter  is  the 
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blast  wave  generated  in  two -phase  detonations  (Ref  10, 11)  and  shock 
ignition  at  liquid  mono  propellant  droplets  (Ref  12).  A  certain  class 
of  time  dependent  energy  input  is  known  to  be  amenable  to  a  similarity 
safaaice  (Ref  2).  Freeman  (Ref  13)  treats  this  problem  in  some 
detail  but  confines  his  treatment  to  a  cylindrical  geometry,  since 
nc  was  interested  mainly  in  applying  his  results  to  spark  discharges 
(Ref  H).  This  report  will  be  concerned  with  the  variable  energy 
blast  wave  and  its  application  to  two -phase  detonations. 
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BLAST  WAVE  EQUATIONS 


The  conservation  equations  applicable  to  blast  waves  are  as 
follows: 


+  IP 

/  3  u 
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for  mass,  momentum. and  energy,  respectively.  Here  transport 
terms  as  might  appear  in  the  momentum  and  energy  equations  are 
omitted.  The  symbol  a  characterizes  the  geometry  considered. 
Thus, 


3  =  0  — *  planar  case 

3=1—*  cylindrical  case 


a=2->  spherical  case 
Combination  of  equations  (3)  and  (1)  yields 


+  u  =  -  yp  f  *iL_  +  ajA  (4) 

dt  ~br  'dr  r  / 

If  the  position  of  the  shock  at  time  t  after  the  initiation  of  energy 
input  is  R  from  the  input  location,  then  the  velocity  of  the  shock  is 
given  by 


U  = 


dR 

dt 


(5) 
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The  boundary  conditions  that  can  be  used  to  solve  equations  (1),  (2), 
and  (4)  are  the  normal  shock  conditions.  Thus,  since  the  velocity 
of  sound  in  the  undisturbed  medium  is 


=  hpoho 

and  the  Mach  number  of  the  shock  wave  is 


(6) 


M  =  U/aQ 


(7) 


the  boundary  conditions  at  r  =  R  are 


(8) 


JL  = 
Po 


Y-l 

Y+l 
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?o 


(Y-l)M2 
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(10) 


These  conditions,  together  with  the  conditions  relating  the  energy 
input  to  the  kinetic  energy  and  the  internal  energy  change  in  the  gas 
behind  the  shock,  are,  in  principle,  sufficient  for  a  complete 
solution  of  the  problem.  The  latter  condition  can  be  written  as: 
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2 

pu^  + 


p-pA  radr 
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(11) 


where: 


Eq  =  Energy  input  per  unit  area 


— .  (Energy  input  per  unit  length) 
2tt 


(Energy  input) 


The  blast  wave  equations  can  be  transformed  in  terms  of 
two  independent  variables  which  are  chosen  to  fit  the  similarity 
solution.  These  variables  are 
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X 

(12) 

and  y,l/M2 

(13) 

The  dependent  variables  then  are 

u  =  U  f(x,y) 

(14) 

P  =  g  (x,y) 

(15) 

y 

:  =  Po  hix,y.» 

(16) 

where  f,  g,  and  h  are  to  be  determined.  If  a  decay  coefficient  is 
defined  as 


x  = 

d.(.nR 

2 

dr_nR 

9  Rd2R/dt2 
TdR/dt)r 


(17) 


Then  equations  (1),  (2),  and  (4),  after  using  equations  (12-17),  can 
be  transformed,  respectively,  to 
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The  boundary  conditions  of  equations  (8)  to  (10)  are  changed  to 


f(l,y)  = 


Y+l 


y .  >  2y 

gd.y)  =  — 


(l-y) 
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Y+l  v+l 


h(l 
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When  an  explosion  scale  radius  is  defined  as 

1 

R0  5  (Vp0>  a71- 

then  condition  (11)  becomes 


(21) 

(22) 

(23) 

(24) 


y(R  /R)a+i  =  J  - 


(a+1)  (Y-l) 


(25) 


where 


( I  “2  +  -yrr)  xadx 


(26) 


It  can  be  seen  that  when  »  1,  y  — >  0  and  therefore  equations 
(18-20)  become  ordinary  differential  equations.  Furthermore,  the 
boundary  conditions  (21-23)  are  simplified,  and  the  last  term  in 
equation  (25)  can  be  neglected.  Despite  these  simplifications, 
however,  the  integral  J  must  still  be  numerically  calculated  after 
numerical  calculation  of  f,  g,  and  h.  The  zero  order  value  of  J, 

i.e. ,  JQ,  under  these  simplifications,  is  found  to  be,  for  Y  =  1.4 
(see  Sakurai,  Ref  2) 


J  =0.596 
o 

(spherical) 

=  0.878 

(cylindrical) 

=  1.696 

(planar) 
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(27) 


It  should  be  mentioned  that  these  values  are  applicable  for  the  case 
where  E„  is  instantaneous.  Under  this  condition,  equations  (25)  and 
(24)  show  that 


R^1  = 


(28) 


and  since  the  independent  variable  y  from  equation  (13)  is 


(29) 


then  by  combining  (28)  and  (29)  and  integrating. 


R 


E 


a  o 

PoJo 


it  is  found  that 
1 

5TT 


(30) 


and 


Eaao 

Vo 


1 

a+3 


(31) 


These  are  the  functional  relationships  of  the  shock  position  and  its 
velocity  for  the  zero  order  solution.  The  decay  coefficient  \Q  in 
equation  (17)  can  be  found  to  be 


X  =  a  +  1 
o 


(32) 


The  decay  coefficient  for  y  >  0  has  been  calculated  by  Sakurai 
(Ref  2),  Freeman  (Ref  13)  and  Goldstine  and  von  Neumann  (Ref  15). 
In  this  case  the  flow  field  does  not  exhibit  similarity,  and  resort 
to  first  and  higher  order  solutions  must  be  made  for  its  accurate 
description  (Ref  2,  13). 
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VARIABLE  ENERGY  INPUT 


As  was  mentioned  before.  Freeman  (Ref  13)  analyzed  in  detail 
the  case  for  a  constant  rate  of  energy  input  in  a  cylindrical  geometry. 
In  this  section  the  zero  order  solution  for  the  general  case  will  be 
presented  and  compared  with  the  instantaneous  energy  input  case. 

For  this,  the  energy  input  is  assumed  to  vary  with  time  as  follows: 


E  =  W  tB 
a  a 


(33) 


where  Wa  is  a  dimensional  constant  of  proportionality  whose  dimensions 
depend  on  the  power  6  .  The  case  of  the  instantaneous  energy  input 
is  covered  by  setting  3=0.  If  an  explosion  scale  radius  is  defined 


(34) 


it  can  be  shown  that  the  shock  radius  can  be  written  as 


R  =  Ktn 

(35) 

where 

1 

/W  a2  Na+S" 

K  =(  ) 

\  p  n2J  ' 
o  o 

(36) 

and 

e+2 

n  = 

(37) 

a+3 

Also  the  decay  coefficient  becomes 


_  Z(l-KX-g) 
2+6 
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It  is  to  be  noted  that  when  5=0,  equation  (35)  reduces  to  equation 
(30)  since  Wa  =  and  equation  (38)  reduces  to  (32).  While  JQ 
appears  in  the  general  case,  it  must  be  mentioned  that  its  value 
depends  on  the  flow  field  (see  equation  26).  Since  the  flow  field 
depends  on  3  via  >.  in  equations  (19)  and  (20),  the  value  of  J  for 
constant  a  has  a  different  numerical  value  for  different  6 .  For 
example,  if  v  =  1.4  and  a  =  1  (cylindrical  case),  JQ  =  .378  for  8=0, 
and  JQ  =  .561  if  8  =  1.  Freeman  (Ref  13)  calculated  JQ  fcr  a  =  1, 
6=0  and  1  for  various  * ;  but  apparently  JQ  for  the  spherical 
case  has  not  been  calculated  for  3=0  (Ref  16). 

To  see  the  effect  of  energy  input  duration*  we  focus  our 
attention  to  cases  where  the  total  energy  is  constant.  Thus, 
equation  (33)  can  be  written  as 

Es  =  Watf  (39) 


where  t^  is  the  duration  of  input  energy.  Substituting  equation  (39) 
in  (36)  and  using  (35),  we  find 


1  3+2 

a  +3  <i  +  3 

.t 


(40) 


Thus,  for  the  same  6»  as  tj  increases,  for  tne  same  time  sifter  the 
start  of  energy  deposition,  the  shock  radius  is  decreased.  A. 
similar  argument  holds  for  the  shock  velocity.  To  see  the  effect  of 
8,  equation  (40)  can  be  written  as 

e+2 

R  =  (t/t  ifTy  (41) 

where 

R  =  R/ 


Equation  (41)  is  plotted  in  Fig  1  for  i  =  Z  and  several  3  (3?0).  The 
nondimensional  radius  R  decreases  with  incrasing  Srand  because  of 
the  change  in  the  value  of  n  and  the  expected  change  in  JQ,  the 
actual  shock  radius  is  expected  to  decrease  with  increasing  g  also. 

While  we  have  focused  our  attention  on  the  zero  order  solution, 
it  must  be  mentioned  that  techniques  for  higher  order  solutions 
also  are  available.  However,  these  are  somewhat  involved;  but 
they  certainly  should  be  used  in  the  future  when  a  more  accurate 
description  of  the  problem  is  warranted. 


li 


COMPARISON  OF  SOME  EXPERIMENTAL  RESULTS 
WITH  THE  VARLABLE  ENERGY  BLAST  WAVE  EQUATION 


As  was  mentioned  in  ihe  Introduction,  shock  ignition  of  mono¬ 
propellant  droplets  can  result  in  a  blast  wave  emanating  from  the 
wake  of  the  drop.  Figure  2,  reproduced  from  the  work  of  Lu  and 
Slagg  (Ref  12),  shows  clearly  a  distorted  spherical  shock  wave 
surrounding  a  2  mm  ethyl  nitrate  drop  ignited  by  a  Mach  3.  3  shock 
wave  in  oxygen.  It  is  reasoned  that,  in  such  a  phenomenon,  the 
reaction  requires  a  finite  time  for  the  energy  release,  and  therefore 
comparison  of  the  experimental  results  should  be  made  with  the 
variable  energy  blast  wave  solution. 

For  complete  combustion  in  oxygen,  ethyl  nitrate  liberates  345 
Kcal/mole,  which  amount  to  approximately  7  x  108  ergs  for  a  2  mm 
drop.  For  spherical  symmetry,  equation  (40)  can  be  written  as 


_ 

1 

3  2 

R  - 

—  (l\z 

E  *  ./ 

(43) 

%Joj  V 

Using  the  values  of: 

E  =  /  x  10  ergs 

pQ  =  5.  6  x  10 g/cm3  (oxygen  at  Mach  3.  3) 
Y  =  1.4 

Jq  =  .596  (assumed  not  to  vary  with  9) 
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we  obtain 


R  =  1.71(t/l0)2/:> 

for  8=0 

(44) 

=  1.45(t/t.)1/:>(t/l0)^/5 

for  3=1 

(45) 

=  1.29il/ti)2/5(t/l0)2/5 

for  3=2 

(46) 

=  1.18(t/ti)3/:>(t/l0)2/:> 

for  3=3 

(47) 

=  1.  io(t/t.)'i/5(t/iO)2/:> 

for  ?=4 

(48) 

7  /5 

=  1.03(t/t  )(t/lO)“ 

for  5=5 

(49) 

where  R  is  in  cm  and  t,  t.  in  usee. 

i 

Equations  (44),  (47),  and  (49)  were  selected  for  comparison 
with  limited  experimental  results  of  Lu  and  Slagg  (Ref  12) .  In 
Fig  3,  it  can  be  immediately  seen  that  the  experimental  data  cannot 
be  explained  by  the  instantaneous  energy  input  case.  Equation  44, 
and  that  the  best  fit  is  obtained  by  the  case  where  8=5  (Equation  49) 
when  t.  is  taken  to  be  10  usee.  It  should  be  recalled  that  this  equation 
is  applicable  up  to  time  t  =  t*  .  Beyond  that  time,  R  is  expected  to 
drop  off  from  what  the  equation  indicates,  due  to  the  fact  that  the 
energy  deposition  ceases  beyond  that  time.  While  more  data  and 
more  accurate  calculations  of  JQ  are  needed  before  final  conclusions 
can  be  drawn,  the  idea  of  interpreting  the  consequences  of  mono- 
propellant  ignition  by  variable  energy  blast  wave  theory  seems  to 
be  reasonable.  The  behavior  of  energy  release  in  a  chemical 
reaction  is  usually  exponential  in  nature,  and  a  manifestation  of  this, 
judging  from  observations  of  the  behavior  of  the  blast  wave,  appears 
also  to  be  true  for  liquid  ethyl  nitrate. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


It  appears  ihat  analysis  of  the  experimental  data  on  shock 
ignited  liquid  monopropellant  behavior  can  be  profitably  made  on 
the  basis  of  variable  energy  blast  wave.  Extensive  data,  if 
compared  with  more  accurate  theory  than  described  here,  should 
lead  to  estimation  of  energy  release  time  and  hence  to  the  initial 
exothermic  reaction  time. 

It  is  recommended  that: 

1.  The  zero  order  solution  for  variable  energy  blast 
wave  be  obtained  for  various  values  of  3,  by  deter¬ 
mining  JQ. 

2.  The  flow  field  within  the  blast  wave  be  evaluated 
more  accurately  than  obtainable  by  the  zero  order 
solution,  using  first  and  higher  orders  for  J*s  and 

X's. 

3.  The  flow  field  and  shock  velocity  be  assessed  at  times 
after  the  cessation  of  energy  release. 

4.  More  data  on  ethyl  nitrate  and  other  monopropellants 
be  obtained  and  compared  with  the  more  comprehensive 
theory. 
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